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Objectives: 



• to evaluate the scene-to-scene registration perforirfnce of the Landsat-4 
MSS relative to previous MSS systems 


t to evaluate the scene-to-scene registration performance of the Landsat-4 
Thematic Mapper (TM). 


Activity this Quarter: 


• A sunmary report dealing with the use of Landsat-4 MSS data in temporal 
data sets has been printed as NASA/NSTL/ERL internal report #222 (a copy 
enclosed). The report has also been submitted to Photogrammetric 
Engineering and Remote Sensing for possible future publication. 

• Reprocessed 16 September 1983 data set of the New Orleans and Mississippi 
Gulf Coast area received. Initial analysis of the radiometric/geometric 
quality comensed. 


Activity Anticipated Next Quarter : 

• Receipt of second TM data set (perhaps a reprocessed December 1983 scene) 
which would permit scene-to-scene registration work to begin 


Problems Encountered 

• Lack of a second TM data set will preclude scene-to-scene work from taking 
place. 
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ABSTRACT 


The MSS sensor on Landsat-4 Is, In certain performance aspects, 
different from those on Landsats 1 through 3. These differences have 
created some concern In the NASA research community as to whether 
Individual data sets can be registered accurately enough to produce 
acceptable data for multitemporal data analysis. This study examines 
the use of Landsat-4 MSS digital data in temporal data sets, and 
presents a method for estimating temporal registration accuracy based on 
the use of an XtY digitiser and grey tone electrostatic plots. Results 
obtained Indicate that the RMS temporal registration errors were not 
significantly different between the temporal data seta generated using 
Landsat-4 and Landsat-2 data (33.35 meters) and the temporal data set 
constructed from two Landsat-2 data sets (33.61 meters). A derivation 
of the model used to evaluate the temporal registration is Included. 

INTRODUCTION 


The purpose of multitemporal registration, using data from a sensor 
such as the Landsat Multlspectral Scanner (MSS), is to produce a single 
multichannel data set from two or more Individual multichannel data sets 
separated In time. The data from one date are used as a base and all 
other data are fitted or mapped onto them. It Is generally accepted 
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that by producing such a teaporal data set researchers can capitalize on 
changes which occur (or do not occur) over the time frame encompassed » 
thus gaining an added dimension that will lead to Improvements In the 
utility of the data being examined* For example, studies by Badhwar, et 
al. (1982), and Hay, et al. (1982), describe the use of temporal 
agronomic crop signatures which characterize Landsat MSS data values as 
functions of crop phenology or growth state* Such temporal analyses 
have produced results which are significantly better than those obtained 
with single-date data sets* Mergerson (1981) points out that as a 
result of better discrimination, multidate data sets produce improved 
crop area estimates* Representative research In the application of 
multitemporal data for Investigations In Forestry (Lee, 1980), Hydrology 
(Salto, et al* 1982), Ice and Snow Field studies (Dey, 1980), Wetlands 
and Coastal Processes (Lim, et al. 1980), Land Cover Dynamics (Byrne, et 
al* 1980), and Land Use Mapping/Natural Resources Inventories 
(Llchtenegger and Seidel, 1981) have demonstrated the widespread 
Interdisciplinary utility of temporal data. 

Results obtained from the use of temporal data are related to the 
degree to which two or more data sets can be registered one to another; 
the better the registration, the better the expected results ceteris 
paribus . Registration errors can arise from numerous sources. Of 
particular Interest to this study are registration errors which might be 
caused by differences In the sensor being used* 

By design, sensor or platform dlffersnces between the first three 
Landsat satellites were Insignificant and thus MSS data registration 
errors Induced by satellite system differences were minimal* However, 
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Che Landsat-4 MSS la In at leaat two respecta different from ita 
predeceaaora > FI rat, the acan angle la wider (14.9 degrees va* 11.6 
degreea), and aecond, the Inatantaneoua field of view (ground resolvable 
area) is slightly different (82*7 x 57 meters va. 79 x 56 meters)* 
Software at the ground processing facility of the NASA Goddard Space 
Flight Center (GSFC) has been designed Co produce "standardized” 
P-^ormat digital data from Che raw Landsat-4 MSS daCar. However, there 
has been some concern expressed about whether the registration of tem- 
poral data sets, which Include Landsat-4 MSS data, has been adversely 
affected (NASA, 1981). GSFC therefore encouraged research in con- 
structing temporal data from Landsat-4 MSS and Landsat 1 through 3 MSS 
data, through the NASA-sponsored Landsat-D Image Data Quality Analysis 
(LIPQUA) Program, to verify continuity of this feature. As its commit- 
ment to the LIDQUA Program, NASA's Earth Resources Laboratory located at 
the National Space Technology Laboratories (NASA/NSTL/ERL) has conducted 
temporal registration assessment research using Landsat-4 MSS data. 


TECHNICAL APPROACH 
Generation of Temporal I^ta Sets 

The procedures for generating temporal data sets are well documented 
in the literature (see especially Bryant, 1982 for a summary). However, 
a brief discussion of the specific technique as implemented at 
NASA/NSTL/ERL is warranted. 

The first step in producing a temporal data set is to choose a 
"base" and a "map" data set. In the context of this paper, base data is 
used "as is”, l.e., it furnishes the scan line and pixel coordinate 
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syttMB onto which the nep deta will be fitted* The nap data set(a) la 
therefore the set which la contorted (atretched/ahrunk) to overlay the 
baae data aet* After designating the base and map data sets, six to ten 
uniformly distributed points which represent unambiguous scene locations 
are located In both data sets* The scan llne/plxel coordinates of the 
points (for both date sets) are used to compute an Initial set of 
mapping equations. Such equations define the fundamental relationship 
between points located In the map data set a..d identical locations In 
the base data set. 

Based on the initial mapping equations, auto-correlation software Is 
used to locate 100 to 200 additional points which represent common 
geographic locations In both data sets. Initially, points are located 
in the base data set, and a window of data centered -on each point is 
extracted. For each such point in the base data set, the Initial 
mapping equations are used to calculate an expected location in the map 
data set. A sliding window the same sice as used In the base data Is 
then moved around the expected location In the map data set, moving one 
pixel or scan line at a time. Each movement results In the extraction 
of a slightly different digital count data sample and the correlation of 
each map data sample with the base data sample is computed. The pixel 
located at the center of the window (In the map data) with the highest 
correlation to the corresponding window of base data Is retained as the 
most probable location for a matching point. This procedure continues 
until all points located in the base data set have been processed. 

Points located by the auto-correlation software are manually edited, 
since some points may not adequately depict the relationship between the 
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bate and map data aete* Thoae which are found to be questionable are 
deleted, And a final set of mapping equations is generated along %rith an 
RMS registration error that is computed usii^g the retained control 
points. Temporal registration is then accomplished using a piecewise 
linear fit with bilinear resampling. 

Registration Evaluation Methodology 

In order to evaluate the level of temporal registration success, a 
technique has been developed at NASA/NSTL/ERL which relies on the use of 
an X-Y digitizer and electrostatic grey-level plots. lAimvrous randomly 
located subscene areas (approximately 512 scan lines by 512 pixels In 
size) are located in the newly constructed temporal data set. Three 
digital count-to-grey level plot'' are generated for each , subscene area 
chosen, one plot for each of two channels of the original base data, and 
one plot for a single channel of the mapped data set after the map data 
set has been resampled. The channels selected are typically MSS channel 
5 (0.7 to 0.8 micron wavelength region) and MSS channel 7 (0.8 to 1.1 
micron wavelength region) for the base data, and MSS channel 5 for the 
mapped data. Since the plots are generated from the temporal data set, 
scan line/plxel coordinates are identical for all plots of a specific 
subscene area, i.e., the scan llnes/plxels of the mapped data are now In 
terms of base data coordinates, since the remapping has already taken 
place. 

For each suhscene area used, the corresponding plots are 
simultaneously mounted on the table of an X-T digitizer. kEter 
initializing the digitizer to the plots, overlay assessment (OA) points. 
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which are unambiguously Identifiable on all three plots, are located and 
their acan llne/plxel coordinates are recorded using the digitiser* For 
each OA point located, three aets of scan llne/plxel coordinates are 
generated, since three plots are used. All coordinates are stored by 
sample number for subsequent use* 

The digitised coordinates are then processed by statistical analysis 
software, which computes the statistics defining the RMS error 
associated with the temporal data set examined. The software Is 
unique In that It computes total error, errors associated with human 
location of OA points, and an estimate of actual temporal registration 
error (which can be thought of as total error minus the human error). 
The following discussion is an amplification of the analjrsls model 
originally designed by Seyfarth and Cook (19R1). 

The model employed assumes that for each point sampled: 


bll 


“ *i ^bll 


where 1 b pixel coordinate for the 1^^ OA test point in the 
first channel of the base data set used, is the theoretically true 
location of the l**^ test point, and is the human error associated 
with selecting that test point. This assumption applies to the second 
base channel used as well as to the mapped channel and yields: 

*1,21 ■ *1 + *b2f 


% 


X - X, + E + a. 
mi 1 ml i 


,th 


where X^21» *ml * coordinate# for the 1 point from the second 

base channel used and the mapped data channel respectively, 
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QummL F.'-tCsU IS 

OF POOR QUALITY 


represent the hunen error contribution for the 1^^ point selected from 

the second base chsnnel used end the napped date channel respectively, 

and Oj is the actual nisrepist ration at the 1^^ point. It is also 

assuned that F^2> independent random variables, 

2 2 

normally distributed with means of tero and variances o o 


, and *nd that a* 


I, and o are all equal to o' 


xm' a' xbl* '' xb2' " xm “ xh 

which is associated with the human error component. (All through this 

discussion similar formulae could be derived for the scan line (Y) 

component, but for the sake of brevity such derivations will not he 

explicitly made.) 

Since and represent coordinates of the same point chosen 

from two channels of the base data set, consider: 


2 ^\li “ *b21^ 

1-1 n 

which can be expanded to yield: 


% 


2 "l * 'bll - *1 - 

i-1 n 

The assumption of independence leads to: 



=b2l2 


® xbl 


- 2 o' 


xh 


2 

xb2 
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ORIGINAL tif 
OF POOR QUALITY 


So (m eetiiMte for the huiMn error coaponent ie: 

• i 2 ‘\u - *b21>' 

1-1 n 


The coordinate derived from the map channel 1* ueed aa followe* 

Conalder the difference between the and the average of X^^^ and 

*b2l’ 


E I 

i-1 


- 


^^11 




■ ^ E pt * *-i * ‘ 'Wi * * 'm‘’1 


i-1 1. 


i £ 

1-1 L 


” ^^'bll ■ *^b21 


21^ 
2” . 


The asaumption of independence yielda: 


n 2 n 2 


" 2 ” 2 
D "1 * k E 1 


i E 


i-1 n 


1-1 n 


i-1 n 


1-1 n 


which la identical to: 

2 2 2 

» ,h * % ♦ i ’ xh 


After replacing o vlth its earlier derived equivalent and rearranging 
tetma we get: 
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ORIGINAL r-cJL ’ m 
Of PCK*R quality 

*al “ ^*bll ' - I H 

i"l 2 1«1 n 

n 

In order to obtain the error in teraa of metexe (aii^e amt are In 
term of plxela) 

^ • 57>J^ for P-fomat Landaat MSS data* 

X me t e rs ^ 

Aa previously oentloned, a alnllar formula could be derived for 

* • Thus, the total RMS error (In meters) would be: 

y meters - 

0 * t ■ V + 0^ 

model ' X meters y meters 

It Is thus possible to obtain estimates of the actual misregistration 
error In temporal data sets. 

Data Sources 

Four P-Format MSS data seta were used in this study as defined In 
Table 1* All scenes Included the New Orleans, Louisiana, metropolitan 
area (Landsat-A worldwide reference system Path 22, Row 39) as well as 
portions of che Mississippi Gulf Coast. A diversity of land covers/land 
use was present including cropland/pasture, dryland, forest, forested 
wetlands, manh, vater bodies, small towns, rivers, etc* 


tandsat 

MSS 

Date 

Scene ID 

Comments 

2 

Sept 

1980 

22063-15500 

Served as the base data set for the study 

2 

Feb 

1981 

22225-15485 


4 

Sept 

1982 

40062-15591 

Scene produced using pre-launch aoftware 
parameters at NASA/Goddard Space Flight 
Center 

4 

Mar 

1983 

40238-16010 

Modified software parameters used 


Table 1. Psta Smirces for the Temporal Registration Investigation. 
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Landsat‘-2 MSS data for 1980 served as the base data set throughout 
study, as they were already In house. All other data were used as 
nap data sets, and were Independently registered to the 1980 data set. 

ANALYSIS ANP RFSIJLTS 

Actual analysis began with the temporal registration of the 1981 
Landsat-2 data set to the 1980 Landsat-2 data set. The results obtained 
(Table 2) were typical of temporal registrations obtained in the past at 
HASA/NSTL/ERL using Landsat 1, 2, or 3 data sets. Subsequent analysis 
of the errors on an OA polnt-by-OA point basis showed no trends in the 
registration errors in tem» of direction of error or magnitude. Thus, 
a particularly "uniform" fit was achieved. 


COMPONENT 

FRROR 

LANDSAT-2 

(1981) 

LANDSAT-4 
' (1982) 

LANDSAT-4 

(1983) 

Pixel Error • 

20.42 meters 

40.32 meters 

20.04 meters 

Scan line Error 

26.70 meters 

26.95 meters 

26.65 meters 

Model RMS Error 

33.61 meters(74)* 

48.50 meters(65)* 

33.35 meters(72)* 


Table 2. Results Obtained for the Registration of Three Scenes of 
Landsat P-Format MSS Pata to a Common Base Data Set. 


*Values enclosed in parentheses represent the number of points used 
to generate the corresponding statistic. 

The 1982 Landsat-A MSS data set was received In January 1983, and an 
Independent temporal registration was produced. The results (Table 2) 
indicated a substantial pixel component error (almost twice the error 
obtained in the Landsat-2 to Landsat-2 registration). The model error 
(estimate of the actual misregistration) also demonstrated a 
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considerable degenerationt undoubtedly the result of the large pixel 
error. When exanlned on an OA polnt-by-OA point basis however, a 
sinusoidal pattern was observed when the magnitude of the pixel errors 
was plotted versus pixel location. Since the position of the MSS scan 
mirror is directly related to the pixel location, this result suggested 
that the coefficients used to model the scan mirror motion were not 
correct. No detectable pattern was noted for scan line errors. 

Tn late January 1983, just before the MSS system onboard Landsat-4 
was turned over to NOAA, corrections were made to the scan mirror 
coefficients used in the Goddard Space Flight Center's ground processing 
software associated with the Landsat-*4 MSS system. The Impact of this 
change can be seen in Table 2. The 1983 data set was collected and 
processed after the new coefficients were incorporated into the 
software. As can be seen, the pixel component error has been reduced 
considerably and is now at the same level as the Landsat-2 only temporal 
registration. e s an line component seems to have been unaffected, 
and remains nearly constant for all of the temporal registrations 
produced, showing no discernible trends in either magnitude or 
direction. 

In order to determine the significance of the changes noted, F-tests 
were conducted on the model RMS errors. The results Indicated that at 
the 10 percent level of confidence, the temporal registrar, on using 1982 
Landsat-4 MSS data (erroneous scan mirror coefficients) was signifi- 
cantly different frcnn each of the two other temporal registrations. The 
results based on the use of the 1983 Landsat-4 MSS data (corrected 
coefficients) were not significantly different from those obtained using 
the 1981 Landsat-2 MSS data. 
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CONCLUSIONS 


It has been denonstrated that temporal registration of Landsat-4 MSS 
P-format digital data with historical data obtained by Landsat-2 can be 
performed to a degree of precision equal to that obtained when using 
Landsat-2 data as both the base and map data sets. Landsat-4 MSS data 
collected between July 1982 and Feburary 1983 should be used with 
caution, since scan mirror motion coefficients were not optimum and use 
of the data may produce less than desirable results. 
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